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ABSTRACT 



This thesis discusses the theory behind the gas dynamic laser. 

A computer program that solves the equations is included along with 
typical outputs. A description of a test apparatus is also presented 
with proof of laser action. 
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I. INTRODUCTION 



The operation of any type of laser system depends on the pro- 
duction of a population inversion. Descriptions of the phenomena 
that take place in a laser may be found in [1] and [2]. There are 
several methods of producing a population inversion, one of which is 
by means of a gas dynamic process. The gas dynamic method usually 
incorporates a combination of rapid heating and rapid cooling of a 
gas . 

In the rapid cooling phase, a gas that is initially at a high 
temperature and pressure is exapnded by flow through a supersonic 
nozzle, which causes a sudden decrease in the temperature and pressure. 
The rapid drop in temperature is accompanied by a decrease in the 
population of the upper energy levels, each of which decays with a 
characteristic time constant. In certain gas mixtures the time 
constants are such that a lower energy level depopulates much more 
rapidly than a higher energy level. The result is a population 
inversion of these two states. The method has been successful in 
producing a laser output as shown by Kuehn and Monson [3]. In their 
experiment a large gas reservoir was heated and pressurized. The gas 
was then passed through a nozzle into an evacuated dump tank. Mirrors 
were placed in the sides of the nozzle to form the optical cavity. 

It was the purpose of this thesis to build and test a gas dynamic 
laser utilizing a shock tube assembly to generate the initial high 
temperature conditions of reservoir gas, and to simulate the device on 
the digital computer according to the theory developed by Basov [4, 3]. 
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II. THEORETICAL DISCUSSION 



The fact that a gas dynamic laser is theoretically possible was 
established by Basov, et al. [4]. This work presented the require- 
ments that must be met by a set of energy- level relaxation times if an 
inversion is to be established. The cases of rapid heating and rapid 
cooling were discussed. The ground states of a CO 2 -N 2 mixture were 
found to satisfy the general requirements. Basov, et al. [5] examined 
this system in detail and offered a set of equations to describe the 
phenomenon. Part of the results of their work is a system of equations 
that describe vibrational relaxation. These equations are 
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This gives a population inversion of 
Nr^o, , - N 



AN "[OCTl] [10 0] , , rt\ , s , , v -2 , . . „ . -2 

iT N (cone) (gj (e 3 - ( e]l + 1) (e' 2 + 1) 



where 



cone = concentration of CCL 



z -k 



In the equations, nitrogen is the "a" gas and carbon dioxide is the 
"(3 M gas. The Z _ represent the number of collisions between a specific 
molecule of type i and any j type molecule. is the probability of 

resonant collision between a CC^ and an Ng molecule. is the corres- 
ponding function for two CC^ molecules. ^ are the probabilities of 

transition from one quantum state to another. The n and m correspond 
to the different energy levels shown in Figure 1. A listing of the 
(n,m) that goes with a given transition may be found in [5]. The 

are the characteristic temperatures of the corresponding levels, and 
hv. 

are given by 0^ = — — . Here, hv^ is the energy of the given state 
above the ground state. The above relations give the rate of change 
of energy in specific energy levels of CC^ and . E^ through 

represent the energies stored in the specific energy level. The initial 

0 /T 

values of E^ through E^ are set by the relation E^ = (Rfh/PL )/(e l -1) 

[5, 6]. 



A basic explanation of the origin of these equations starts by 
considering the change of population of some state of one particle as 
the particle undergoes collisions with another particle. The rate of 
change of population would be given by 
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where one is now considering an N- level system. P is the probability 

that a particle will undergo a transition from the n^ energy level to 

the m^ energy level. (Basov uses p^ to indicate a probability of 

transition due to resonant collisions and to indicate a probability 

due to other than resonant collisions.) The first term on the right- 

hand side gives the depopulation of the i t ^ 1 level. The second term 

introduces the re population of the i level. Z is the number of 

e,p 

collisions that an e type particle will undergo with a p type particle. 
Note that only transitions due to interaction with another particle 
are being considered at this point. Transitions that may be caused 
by other reasons (spontaneous, photon stimulated, etc.) are not included. 
These were assumed negligible. 

In order to describe the population inversion process, one must 

know the transition probabilities (p ) and the collision rates Z 

r r nm e,p 

The actual calculations that must be carried through to obtain the p 

r nm 

will not be worked out here. Rather, an intuitive feel for the method 
used will be presented. One starts by writing Schroedinger ' s equation 
for the collision process [6]. This will result in a wave- function 
description of a particle that is under the influence of an intermole- 
cular interaction potential. This interaction potential is a function 
of the distance between the involved particles. 

The resulting wave functions can be written down as some combination 
of the wave functions of an isolated particle. The set of constants 
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that perform this transformation is directly related to the transition 

probabilities p . The relationship between the elements of the set 
nm 

of constants and the transition probabilities is known, and is of 
sufficient simplicity as to allow calculation of the probabilities. 

The origin of the *s is discussed in [5] where they are derived 

from the kinetic theory of rigid molecules. The probabilities given 
by Basov are altered versions of the theoretical relationships so that 
the best experimental correlation results. 

In order to form a closed system of equations that describes the 
gas dynamic process, one must describe the temperature, pressure and 
density as functions of time. This is necessary due to the dependence 
of the probability relationships on these parameters. The quasi one- 
dimensional flow equations form a set of independent equations that 
describe the above variables. One of the relationships is the equation 
of state: 




where R is the universal gas constant and M is the molecular weight of 
the gas. If there is more than one gas, this equation can be separated 
into the partial pressure components. That is, 



Assuming steady-state flow, 
equation, which is given by 



P. = D . r- T. 
1 1 M. 

l 



another relationship 



is 



the continuity 



^ (pav) =0. (6) 

A third relationship is that for the conservation of momentum, as given 
by 
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dx 
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Finally, there is the conservation of energy for adiabatic flow. This 
equation takes the form: 






( 8 ) 



and in effect states that there is only a finite amount of energy that 
can be distributed among the different modes of energy storage. 

These eight equations allow one to describe the relaxation pehnomenon 
as exhibited by an adiabatica lly cooled gas. A digital computer program 
that solves the equations is shown in the Appendix. 

Before describing the theory of the gas dynamic process, keeping 
the theory of the ^-CC^ gas dynamic laser in mind, the ^-CO^ mole- 
cular system will be discussed. The 10.6-micron molecular laser 

system inverts the population of certain of the molecular vibrational 
energy levels of the ground state of the CO^ molecule. The vibrational 
energy levels are split by the molecular rotational energy. However, 
in the energy-level diagram, for clarity purposes, these rotational 
energy levels are not shown. The different vibrational modes of the 
ground state of the CO^ and molecules are shown in Figure 1. 

With reference to Figure 1, is the lowest vibrational level of 
the atom. Since ^ is diatomic, there is only one vibrational degree 
of freedom. This is a symmetrical stretching mode that can be visualized 
as similar to two balls that are attached by a spring. 

(00°1) is the energy of the assymetrical stretching mode of the CO^ 
molecule. This is the case of the two oxygen atoms moving in one 
direction while the carbon atom moves in the other" direction . 

(01 0) is the bending mode energy. This is a degenerate state, as 
the molecule can bend in either of two perpendicular planes, the degeneracy 
being indicated by the superscript. 
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(10°0) is the energy of the symmetric stretching mode. This is 
the same as in the case of ^ , except that there is a "stationary” 
carbon atom between the vibrating oxygen atoms. 

The above energy levels are simple harmonic oscillator states, and 
therefore the energy difference between the nth and (nth+1) energy 
levels is the same, regardless of the value of n. The nomenclature 
for the higher vibrational energy levels is the same as above. For 
example, the level (02°0) is the nondegenerate bending mode that has 
an energy of approximately twice the (01^0) level. Note that although 
it is not shown, the (02°0) would be right below the (10°0) level. 

There are other levels in the region shown by Figure 1, but the 
levels shown are those that enter into the strong 10.6-micron transition, 
and are the most important ones for the present purposes. 

With the energy levels set forth, it is now necessary to explain 
several other phenomenon before the laser can be completely described. 
The first of these phenomenon is what is known as a collision of 
the second kind or a resonant collision. Basically the resonant 
collision is the interaction of two atoms, one of which is in an 
excited state. The other atom may be in any state. The result of 
the interaction is that the excited atom transfers its energy to the 
second atom, while returning itself to the ground state. This may be 
written schematically as [2]: 

X + Y* => X* + Y 

where the * indicates an atom in an excited state. It has been shown 
that the probability that the interaction will occur decreases as the 
energy gap between the energy level of atom Y and the closest available 
energy level of atom X increases. If the energy difference is greater 
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than 800 cm , the collision will effectively not occur [7]. Although 
a proof will not be presented here, some insight can be gained by 
considering the fact that both the energy and the momentum of the two 
atoms must be conserved. This criterion is most readily fulfilled when 
there is little "internal" energy (i.e., electronic, vibrational, or 
rotational) that must be transformed into kinetic energy of trans- 
lation of the atoms. 

The next item that must be discussed is the ^ energy- level diagram. 
As seen from Figure 1, for ^ there is only one state of interest: the 
lowest vibrational level of the ground state. This state is a molecular 
level, and since has no net electric dipole moment, there can be no 
transitions directly to the ground state due to any electromagnetic 
interaction. The result is that for the ^ molecule to return to the 
ground state from the vibrational levels it must undergo a collision 
with another particle (other molecules or the container walls). This 
is known as a metastable state, and is characterized by a long lifetime. 

With the above information, the CC^-^ molecular gas dynamic laser 
system can be described. Once the gas has been heated, populating the 
upper levels, it is then expanded by means of a supersonic nozzle. 

The excited ^ molecules interact with the CO^ molecules, and the two 
types of molecules undergo resonant collisions which serves to maintain 
the population of the 00°1 level of the . The population of the 
other CO 2 energy levels depopulate rapidly as both levels (shown in 
Figure 1) are connected to the ground state by large transition 
probabilities. The large energy gap between these states and pro- 
hibits repopulation due to resonant collisions. This creates a 
population inversion between the 00°1 and the 10°0 states. There are 
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other inversions, for example between the 00°1 and 02°0 levels, but 
the above inversion is the strongest, and gives the 10.6-micron 
transition. The lasing action may now occur. 

There is one other characteristic of the laser system that 

should be pointed out. This is the fact that any impurities in the 
gases will cause a marked decrease in the output power. This is due 
to the resonant transfer of energy from the ^ molecules to the 
impurity atoms, which robs energy that would have otherwise been 
pumped into the inversion. 

Although the above theory does describe the same ^-CC^ laser 
system, higher output powers can be obtained by the use of additives. 
In the case of a He-^-CC^ laser, a major factor in the population of 
the 00°1 level appears to be due to cascading effects from higher 
energy levels of both molecules [8, 12]. That is, there is resonant 
transfer of energy between the CC^ and He at higher energy levels 
(collisions of the second kind at upper levels). The CC^ molecules 
then spontaneously drop into the 00°1 level to set up the inversion. 
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III. EXPERIMENTAL APPARATUS 



The previous section dealt with the theory behind the gas dynamic 
laser. The apparatus used to test the theory is shown in Figures 2 
through 13. The shock-tube facility used had been previously built 
at the Naval Postgraduate School and it is discussed in detail by 
Penaranda [13]. 

Figure 2 shows a schematic of a double-diaphragm shock tube. The 
shock tube generates a shock wave by suddenly expanding a high-pressure 
gas into a low-pressure gas. The shock wave may be visualized as 
formed by the superposition of inf initess imal compression waves. The 
shock wave is accompanied by large gradients in temperature and pres- 
sure. The double-diaphragm arrangement was used because of the repro- 
ducability of run conditions that it affords. In the double-diaphragm 
assembly, one drops half of the driver pressure across each diaphragm. 
The diaphragms (flat brass plates) are cut to withstand this pressure. 
However, when valve #1 is opened the entire driver pressure is felt 
across diaphragm #1 which has been cut to rupture well below this 
pressure. Therefore, the shock will always be originated by the same 
total driver pressure. This system eliminates the necessity that the 
diaphragms rupture at an exact pressure. The only requirement is that 
they rupture in the pressure range between the full driver pressure 
and one half that value. 

Figure 3 shows the regions of interest for the case of a reflected 
shock wave. The shock wave is traveling into the driven gas (region 1) 
at a velocity S. The contact surface is propagating into the driven 
gas at a velocity C. This surface separates Region 2 from Region 3. 
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Regions 1 and 2 contain pure driven gas while Region 3 contains pure 
driver gas. That is, there is no mass flow across the contact surface. 
The array of lines in Region 4 represents the rarefaction wave, which 
is produced by the expansion of the driver gas. Region 4 is the driver 
gas at rest. Region 5 contains pure driven gas after the shock has 
reflected from the end wall of the shock tube. The conditions of the 
gas in Region 5 were used as the initial conditions for the experiment. 

Analytical descriptions of the state of the gas in any of the 
regions may be obtained by means of the conservation of mass, momentum 
and energy equations in conjunction with the adiabatic gas laws. This 
is done in (9) and (10). The results of interest to the present experi- 
ment are 



' 27 M 2 - (7-1)1 [ (37-1) M 2 - 2 

(7+1) J L ^ ^ m 2 + 2 

^5 = [2 (7-1) M 2 + (3-7)1 f (3-7-1) M 
T 1 [(7+1) M] 2 





- 2 



and 



(7-DI ( 



where y is cp/cv, a^ is the speed of sound in Region 1 and M is the 
mach number of the shock wave. The conditions in Region 5 can be 
controlled by variation of M, P, and a^/a^. 

The remainder of the physical apparatus is shown in Figures 4 
through 12. Figure 4 is a schematic of the entire setup. The three 
amplifiers shown serve to amplify the signals from the thin film 
gauges (T.F.G. in the diagram). The distance between T.F.G. #1 and 
#2 is 5 ft . 3 in. Therefore by recording when the shock passes each 
of the gauges, the actual shock speed can be determined. The schematic 
of the amplifiers is shown in Figure 5. 

The biasing circuit for the detector is shown in Figure 6. The 
detector is biased at about 370 microamps and gives a signal of several 
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millivolts. The detector is gold-doped germanium with a bandpass of 
1.0 to 10.6 microns. The time constant of the detector is 50 
nanoseconds . 

The ’’dump tank" serves as a reservoir for the gas after it has 
passed through the optical cavity formed by the mirrors. As will be 
seen in a latter portion of this section, the observation time is 
on the order of 2 milliseconds. The particle velocity will be super- 
sonic after the gas expands through the nozzle. For a particle 
velocity of 4000 fps , the dump tank must be 4 feet long to avoid 
interference of the returning gas with that gas that is still expand- 
ing through the nozzle. The gas that reflects from the end wall of 
the dump tank will introduce a change in the conditions within the 
gain path. This change in the conditions marks the end of the test 
time. The dump tank can be seen in Figure 7. The fingerlike pro- 
jection from the endwall is a safety release valve, which is set to 
release when the pressure exceeds 200 psi. This is well below the 
rated pressures of all components in the system and will therefore 
prevent any overpressures from building up. 

The mirror and nozzle assembly can be seen separated in Figures 8 
and 9. A mirror mount is shown disassembled in Figure 10. A cutaway 
view is shown in Figure 11. The head of the bolt (shown in Figure 10 
with the mirror attached) is ground off so that the mirror surface is 
recessed 0.075 inch from the nozzle wall. This is to help preserve 
the mirrored surface. Original test runs were made with gold-coated 
glass mirrors (not recessed) and mirror surface deterioration was 
severe. In some places the gold was dislodged from the glass. The 
recessing of the mirrors alleviated this problem. 
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The mirror has a radius of curvature of 4.5 inch (2.25-inch focal 
length). The gain path with the mirrors recessed is 3.14 inches. The 
mirrors are 98$ reflective at 10.6 microns. 

The nozzle is a linear nozzle with a half angle of 30°. The throat 
is 2 millimeters high. The center of the mirrors are 1.0 inch down the 
nozzle from the throat. The inside of the nozzle is roughened in order 
to minimize internal reflections of spontaneously emitted photons. 

The plexiglass shield shown in Figure 12 absorbs any stray 
radiation. The end wall is at 10° to the incident beam (which enters 
from the left of Figure 12). This will relfect 7$ of the energy. The 
reflecting surface is an IRTRAN 2 crystal. The hole drilled in the 
plexiglass behind the IRTRAN is to isolate the IRTRAN from the load 
(that is from another piece of plexiglass that will absorb the trans- 
mitted beam). The two metallic bars on the far side serve as counter- 
balances to the detector (which is filled with liquid N^ under operation). 
The white wire carries the signal while the black lead grounds the 
detector and Dewar flask. 

Figure 13(a) shows the output of the thin film gauges of the shock 
tube for a test run (0.5 m sec/cm). The conditions were 600 psi N^ 
driving 40 mm Hg. ARGON was the driven gas. The resulting mach speed 
was 4.0. An endwall pressure gauge was substituted for the nozzle and 
its output is shown in Figure 13(b) (0.5 m sec/cm). The flat region in 
the middle of the trace is the observation time for Region 5, and is 
found to be on the order of 2 m sec. 



24 



IV. RESULTS 



A. THEORETICAL 

Figures 14 through 18 show different outputs from the computer 
program. Figure 14 shows a comparison of the program solutions to 
the predictions by Basov, et al. [5]. The case considered was for 
an axially symmetric nozzle, with initial test conditions of 10 
atmospheres and 1500° Kelvin. The two solutions were found to differ 
by no more than 3$. 

Figure 15 shows a comparison of the linear nozzle as used in the 
experiment to the hyperbolic nozzle considered by Basov. It should 
be noted that the pressure drops faster in the hyperbolic nozzle. 

The effect that this has on the population inversion can be seen in 
Figure 16. The two larger curves show comparisons of the computer 
predictions to Basov's predictions for the hyperbolic nozzle. The 
smaller population inversion in the linear nozzle is caused by the 
higher pressure and temperature (as compared to the hyperbolic nozzle) 
Also, since the temperature decreases more slowly in the linear 
nozzle, the peak in the population inversion occurs farther down the 
nozz le . 

Figure 17 shows the population inversion at the center of the 
optical cavity as a function of temperature. There is some optimum 
temperature for each initial pressure and CO^ concentration. Below 
this temperature less energy is pumped into the upper levels. Higher 
temperatures increase the gas velocity through the nozzle. Since the 
characteristic relaxation times are effectively independent of tempera 
ture in this range, the population inversion peak occurs farther down 
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the nozzle than the center of the optical cavity. This causes a 
smaller inversion at the center of the optical cavity. 

B . EXPERIMENTAL 

Figure 18 shows the result of one of the initial tests. Time 
limitations precluded further experimentation. The setup used in 
this initial test for lasing action was a modification of that shown 
in Figure 4. The modified setup used only the dual-beam oscillo- 
scope. The oscilloscope was triggered by the first thin-film gauge. 

The second and third gauges were fed into different (i.e., opposite 
polarity) differential inputs of one channel. This was the lower 
trace in Figure 18. The detector output was fed directly into the 
second channel. The positive jump in the lower trace at 3.8 msec is 
the difference of the outputs of the second and third thin-film 
gauges, and comes just prior to reflection of the shock wave off the 
endwall. At this time the detector output shows a negative step, 
which was caused by the lasing action. Two milliseconds later the 
conditions of the gas were upset and the observation time as described 
by the computer program was over. This is noted by a slight ripple 
in the upper trace at 6.0 msec. The large spikes are due to recreation 
of favorable conditions by internal shock reflections. Their large 
size is due to mixing of the helium driver gas. This was shown by the 
fact that the spikes disappeared (but not the original step) when 
was used as a driver gas. 
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V. CONCLUSIONS 



From the reproducibility of the output from the detector, it was 
concluded that lasing action was taking place. However, insufficient 
data were taken to permit comparison of the actual laser output to 
the theoretically predicted value. The only indication that the theory 
correlates with the experiment is the fact that a population inversion 
was predicted for the test conditions that resulted in laser output. 

Further work with the system, as shown in Figure 4, is necessary 
to supply the data required for full comparisons between the theory 
and experimental results. Further experimentation with the laser 
should supply insight into both the validity of the theory and into 
the shock reflection and gas mixing that caused the large output 
spikes . 
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APPENDIX 



The eight equations given in. Section I form the closed set 
necessary to describe the gas dynamic laser. The computer program 
at the end of this section solves these equations. Descriptions of 
the variables may be found in tables I and II. There are a few 
portions of the program that need clarification. 

The factor (TT/TC) in statement 25 ensures that T = TT during the 
first iteration. If T ^ TT, a divergent solution results. In some 
cases the solution may converge after several iterations. The error 
between T and TT is small and is due to round-off errors of the 
conversion factors used. 

2 

The factor (0.98) in the relation for CSQ (C ) also ensures 

computational stability. That is, in order to compute dp/dx, one must 
2 2 

divide by (1-v /C ) . At x = 0 (at the throat), v = c, and one would 
therefore be dividing by zero. The factor 0.98 ensures that this 
situation does not occur. The resulting error is less than Vf>. 

The subroutines RKLDEQ and OSPLOT may be found in the USNPGS 
System 360 program library. More information on them may be obtained 
there . 
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TABLE I 



Variable 


Physical Interpretation 


Units 


Y(1),Y(2),Y(3),Y(4) 


Total energy in given level 


cm V gr 


P5, PI, V(5) 


Pressure 


atmospheres 


P1AM 


Actual driven pressure 


millimeters 

mercury 


TT,T1,T5,T,TC 


Temperatures 


degrees Kelvin 


Y (6) , VT 


Gas particle velocity 


cm/sec 


CSQ 


Speed of sound squared 


2, 2 
cm /sec 


F(1),F(2),F(3),F(4) 


Population rates 


cm" V gr/ sec 


F(5) 


Pressure rate of change 


atm/ sec 


F(6) 


Velocity rate of change 


cm/sec/ sec 


PI 


Normalized population 
inversion 


dimens ionless 
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TABLE II 



Computer Variable 


Basov Variable 


Y(I),Y(2),Y(3),Y(4) 


E l’ E 2 ' ’ E 3’ E a 


Y(5) 


P 


Y(6) 


V 


F (i) 


F(i) 


T5 


T 

o 


P5 


P 

o 


CONC 


7 


Z. . 


Z. . 


Q ij 


Q ij 


P. . 


P. . 
iJ 


PI 


AN/N 

o 
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FIGURE I 

C0 2 AND N 2 ENERGY LEVELS (in cm 1 ) 
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FIGURE 2. 

DOUBLE DIAPHRAGM SHOCK TUBE. 
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Shock Tube At Time t, >0 
R C S 




FIGURE 3 

PLOT OF DISTANCE VS. TIME IN SHOCK 
TUBE AFTER DIAPRAGM IS BROKEN. 
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Shock Tube End Wall 



nMirror And Nozzle 
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FIGURE 4 

SHOCK TUBE LASER DIAGRAM 



ADJUST NULL BIAS WITH A VTVM BEFORE 
CONNECTING OUTPUT TO AMPLIFIER. 




FIGURE 6 
BIASING. CIRCUIT 
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FIGURE 7. 

END VIEW OF SHOCK TUBE 
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FIGURE 8 
MIRROR MOUNT 
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FIGURE 9 
NOZZLE ASSEMBLY 
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FIGURE 10 

MIRROR MOUNT (DISASSEMBLED) 
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FIGURE 12 

PLEXIGLAS SHELD AND DETECTOR 
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TFG # 2 



TFG # 3 




FIGURE 13(a) 



OUTPUT OF T.FG. FOR TEST RUN 
(0.5 m. sec./ cm) 




FIGURE 13(b) 

OUTPUT OF ENDWALL PRESSURE GAUGE 
(0.5 m sec. /cm) 
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FIGURE 14 

COMPUTER OUTPUT COMPARISON FOR IDENTICAL NOZZLE. 
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BASOV PREDICTION 
}C0MPUTER RESULTS 



A 




52 



FIGURE 16 

POPULATION INVERSION VS. DISTANCE DOWN NOZZLE 
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FIGURE 17 

POPULATION INVERSION VS. TEMPERATURE AT 
CAVITY CENTER (LINEAR NOZZLE) 



Detector 



T.F.G. 




2 msec/cm 

Lower and upper trace; 10 mv/cm 
Run conditions; 

Driven gas: 200 ram (Hg) of 20,7# C0 2 in N 2 
Driver gas: 600 psi of U0% Argon in Helium 



FIGURE 13 

DETECTOR OUTPUT SHOWING LASER ACTION 
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